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Modern applications of Covariant Density Functional Theory (CDFT) are discussed. 
First wo show a systematic investigation of fission barriers in actinide nuclei within 
constraint relativistic mean field theory allowing for triaxial deformations. In the second 
part we discuss a microscopic theory of quantum phase transitions (QPT) based on the 
relativistic generator coordinate method. 

1. Introduction 

Theoretical approaches based on Density Functional Theory (DFT) remain un- 
doubtedly among the most successful microscopic descriptions of nuclear structur^. 
One of the r nost successful schemes of this type is covariant density functional the- 
ory (CDFTj2El_ jt is based on Lorentz invariance, connecting in a consistent way 
the spin and spatial degrees of freedom in the nucleus. Therefore, it needs only a 
relatively small number of parameters which are adjusted to reproduce a set of bulk 
properties of finite nuclei and it is valid over the entire periodic table. 

Over the years, Relativistic Mean-Field (RMF) models based on the CDFT have 
been successfully applied to describe ground state properties of finite spherical and 
deformed nuclei over the entire nuclear chart. Excited states have been investigated 
in the framework of relativistic RPA and the relativistic Cranking modeM Here 
we restrict ourselves to two more recent applications of this concept. In section 2 
we discuss a systematic investigation of the inner fission barriers in the actinide 
region and in section 3 we show an application of the relativistic Generator Co- 
ordinate Method (GCM). It allows a microscopic description of quantum phase 
transitions in finite nuclei, which have been investigated extensively in recent years 
by phenomenological models. 



2. Triaxial fission barriers in the actinide region. 

Most of the investigations of fission barriers have been carried out in the frame- 
work of non-relativistic modelJ^. A study of the inner fission barrier heights Bj* of 
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even-even nuclei is motivated by the importance of this quantity for several physical 
phenomena. Many heavy nuclei decay by spontaneous fission, and the size of the 
fission barrier is a measure for the stability of a nucleus reflected in the spontaneous 
fission lifetimes of these nuclei^. The probability for the formation of a superheavy 
nucleus in a heavy-ion-fusion reaction is also directly connected to the height of 
its fission barrier. _B|* is a decisive quantity in the competition between neutron 
evaporation and fission of a compound nucleus in the process of its cooling. The 
large sensitivity of the cross section a for the synthesis of the fissioning nuclei on 
the barrier height S^* stresses the need for accurate calculations of this value. For 
example, a change of by 1 MeV changes the calculated survival probability of 
a synthesized nucleus by about one order of magnitude or even more. The popula- 
tion and survival of hyperdeformed states at high spin also depends on the fission 
barriers. In addition, the fission barriers of very neutron-rich nuclei have (among 
other quantities such as Q-values or /3-decay rates) an influence on the r— process 
of stellar nucleosynthesis. 

During the last decade the role of triaxiality in the region of the saddle point 
of fission barriers has been recognized and tested in many non-relativistic frame- 
works. It was f ound that the height of the barrier is reduced when triaxial shapes 
are allowed^^. However, this lowering strongly depends on the proton and neu- 
tron numbers and on the model employed. Investigations of inner fission barri- 
ers with triaxiality included are available within the frameworks of the micro- 
scopic-|-macroscopic method, the extended Thomas-Fermi plus Strutinsky integral, 
and non-relativistic energy density fun ction als based on Skyrme and Gogny forces. 

Covariant density functional theorjEI^ is an approach alternative to the above 
mentioned non-relativistic methods. The progress in the study of the fission barriers 
within CDFT has been slower than in its non-relativistic counterparts. Inner fission 
barriers in several nuclei have been calculated mostly in the axially symmetric rel- 
ativistic mean field (RMF) -I- BCS approach. Many of these investigations use the 
constant gap approximation in the BCS part. However, our recent study of various 
pairing schemes^ used for calculations of fission barriers clearly shows that this 
approximation leads to unphysical results, because the gap changes considerably 
along the fission path. Only relatively few fission barriers have been studied in axi- 
ally symmetric RMF calculations within the BCS approximation using an effective 
density-dependent zero-range force in the pair ing ch annel presenting a much more 
realistic approximation for the pairing channel^^^El 

Unfortunately, axially symmetric calculations cannot be directly compa red w ith 
experimental data since, as has been shown in non-relativistic calculations!^^, the 
lowering of fission barriers due to triaxiality is significant and can reach 3-4 MeV 
in some nuclei. At present, no systematic studies of the effects of triaxial degrees 
of freedom on the height of inner fission barriers are available in covariant density 
functional theory. Here we discuss recent systematic investigationJI^ of the inner 
fission barriers within the triaxial RMF-I-BCS approach, and confront these CDFT 
results with experimental quantities as far as possible. 
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Fig. 1. (Color online) Potential energy surface in ^''"Pu. The energy difference between two neigh- 
boring equipotential lines is equal to 0.5 MeV. The blue dashed line with solid circles shows the 
lowest in energy solution as a function of /32. Further details are given in the text. 



In this investigation we use the parameter set NL3 which is a modern version 
of the famous set NL3 and we employ in the pairing channel a monopole force with 
a set of coupling constants G(N, Z) adjusted carefully to experimental odd-even 
mass differences in this regioill^. In Fig. [1] we show as an example the potential 
energy surface of the nucleus ^^''Pu in the /3-j plane. For axial symmetry we find 
the normal deformed minimum of the ground state at a deformation f3 0.28, a 
maximum at (3 ~ 0.52 and a superdeformed minimum at /3 ^ 0.96. We observe 
that the fission path (the blue dashed curve) bypasses the axial barrier between 
the normal and superdeformed minima. The barrier height is determined by the 
maximum of the energy along this fission path. 

The deformation energy curves for other even-even nuclei in this region obtained 
in these calculations are shown in Fig. [2] Full black curves show axially symmetric 
solutions, while the values of the deformation energy curves along the triaxial fission 
path are given in red. One can see that by allowing for triaxial deformation the 
fission barrier heights are reduced by 1 — 4 MeV as compared with axially symmetric 
solutions. This lowering depends on the proton and neutron numbers. It also brings 
in average the results of the calculations in closer agreement with experimental 
data shown by green solid circles in Fig. [2] These circles display the height of the 
experimental fission barrie/^ at the calculated /3-deformation of the saddle point. 
A typical uncertainty in the experimental values, as suggested by the differences 
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Fig. 2. (Color online) Deformation energy curves of even-even actinide nuclei obtained in the 
RMF+BCS calculations with the NL3* parametrization. Black solid lines display the deformation 
energy curves for the axially symmetric solution, while red solid lines show the deformation energy 
curves along the triaxial part of the fission path. Solid green circles show the experimental values 
for the height of inner fission barrier. 

among various compilations, is of the order of ±0.5 MeV. 

The microscopic origin of the lowering of the barrier due to triaxiality can be 
traced back to the changes of the level density in the vicinity of the Fermi level 
induced by triaxiality. The single-particle level density at the Fermi level is lower for 
triaxial solutions than for axially symmetric solutions. This is especially important 
at the deformation corresponding to the saddle point of the axially symmetric 
solution. Fig.[3]shows the differences between calculated and experimental heights of 
inner fission barriers. The average deviation between theory and experiment is 0.76 
MeV. This is comparable with the results obtained in the macroscopic-|-microscopic 
method. 

3. Microscopic theory of quantum phase transitions. 

Applications of CDFT discussed so far are based on the mean field concept. How- 
ever, this framework has its limitations. It is characterized by symmetry breaking. 
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Fig. 3. (Color online) The difference between experimental and calculated heights of inner fission 
barriers as a function of neutron number N and proton number Z 



and can only provide an approximate description of specific nuclear excitations, 
such as collective rotations in the framework of the Cranking model or vibrations 
in the framework of time-dependent mean field theory, which corresponds in the 
small amplitude limit to RPA or QRPA. To calculate excitation spectra with spe- 
cific quantum numbers and electromagnetic transition rates in individual nuclei, it 
is necessary to extend this framework to include collective correlations related to 
the restoration of broken symmetries and to fluctuations of collective coordinates. 
Collective correlations are sensitive to shell effects, display pronounced variations 
with particle number, and cannot be incorporated in a single Slater determinant. 

On this level an extension of density functional theory beyond mean flcld takes 
into account collective correlations through the restoration of broken symmetries 
and configuration mixing of symmetry-breaking product states. In recent years 
several accurate and efficient models and algorithms have been developed that 
perform the restoration of symmetries broken by the static nuclear mean field, and 
take into account fluctuations around the mean-field minimum. The most effective 
approa ch to configuration mixing calculations is the generator coordinate method 
based on the variational ansatz for the many-body wave function: 



l*a) 



J dqUq)P'p''P^\^{q)) 



(1) 



The set of intrinsic wave functions |$(<z)} is generated by self-consistent constraint 
RMF+BCS calculations. The weight functions fa{q) in the collective wave function 
of Eq. ([T]) are determined from the variation of the energy by the Hill- Wheeler 
equation: 



dq' [H{q, q') ~ E^N{q, q')] U{q') = 0, (2) 

where H and M are the angular-momentum projected GCM kernel matrices of the 
Hamiltonian and the norm, respectively. With the simplifying assumption of axial 
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symmetry, GCM configuration mixing of angular-momentum, and even particle- 
number projected quadrupole-deformed mean-field states, has become a standard 



Much more involved and technically difficult is the description of intrinsic 
quadrupole modes including triaxial deformations. Only very recently fully mi- 
croscopic three-dimensional calculations have been presented for Skyrmel^ and 
Qognyl^ mean field models generated by triaxial quadrupole constraints that are 
projected and mixed by the generator coordinate method. This method is actually 
equivalent to a seven-dimensional GCM calculation, mixing all five degrees of free- 
dom of the quadrupole operator and the gauge angles for protons and neutrons. 
Only very recently 3D angular momentum projection " and configuration niixing 
of angular-momentum projected triaxial relativistic mean-field wave 

functionPQEll 

has been possible. However, triaxial calculations of this type are connected with an 
extreme numerical effort and are, at present, not possible for heavy nuclei. 

Here we report on an axially symmetric application, where the set | <&((?)) is 
generated by imposing constraints on the axial (720 mass quadrupole moments in 
self-consistent RMF-I-BCS calculations. The following calculations are based on 
the relativistic point coupling model PC-Pi^Sl. Pairing correlations are taken into 
account in the simple BCS approximation, using a (5-pairing force. 

Quantum phase transitions from spherical to deformed shapes have recently 
gained considerable interest. However, most of these investigations have been car- 
ried out within simple boson models^. In order to understand the microscopic 
background of such investigations in finite nuclei the relativistic Generator Coordi- 
nate Method has recently been appliecP^. Starting from constrained self-consistent 
mean-field calculations of potential energy curves, the one-dimensional generator co- 
ordinate method (GCM) has been used to perform configuration mixing of angular- 

0/1 

momentum and particle-number projected relativistic wave function^^. It has been 
shown that the microscopic framework based on universal density functionals, ad- 
justed only to nuclear ground-state properties and extended to take into account 
correlations related to symmetry restoration and fluctuations of collective variables, 
is able to describe not only general features of such shape transitions but also the 
unique behavior of excitation spectra and transition rates at the critical points. In 
particular the shape phase transition in the chain of Nd isotopes occurs precisely 
at the isotope ^^^Nd, in agreement with empirical evidence for the realization of 
the X(5) symmetry in the N = 90 rare-earth isotones. As shown in Fig. [4] the cal- 
culated spectrum and the intra-band and inter-band B(E2) transitions for ^^''Nd 
are in excellent a gree ment with data^'' and close to the X(5)-model predictions for 



As discussed before, because of the considerable numerical effort 3D angular 
momentum projected GCM calculations are at present not yet possible for such 
heavy nuclei as the Nd-chain. Therefore this phase transition has been investigated 
recentl}I2S! on the basis of a 5-dimensional Bohr Hamiltonian, whose parameters 




tool in nuclear structure studies with non- relativisticM^Lk^, and relativistic density 




functionalJlSlII!. 
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Fig. 4. The particle-number projected GCM spectrum of ^^"^Nd (left), compared with the and 
the X(5)-symmetry predictions (right) for the excitation energies, and B(E2) values (in Weis- 
skopf units) of the ground-state and the /3-band. The theoretical spectra are normalized to the 
experimental energy of the state 2^, and the X(5) transition strengths are normalized to the 
experimental B(E2;2+ ). 



were derived microscopically form the covariant density functional PC-Fl. The 
resulting spectrum and the electromagnetic transition probabilities are more or less 
the same as in the ID- AMP GCM calculation shown in Fig. H 

4. Summary and outlook. 

Among the microscopic approaches to the nuclear many-body problem, the frame- 
work of nuclear density functional theory provides the most complete description 
of ground-state properties and collective excitations over the whole nuclide chart. 
Here we have seen that the self-consistent relativistic mean-field model based on 
a universal density functional provides an excellent description of the inner fission 
barriers in actinide nuclei with an average deviation between theory and exper- 
iment of 0.76 MeV. For transitional nuclei relativistic models not only describe 
general features of shape transitions but also particular properties of spectra and 
transition rates at the critical point of the QPT. However, to calculate excitation 
spectra and transition probabilities, the self-consistent mean-field approach has to 
be extended to include correlations related to restoration of broken symmetries and 
fluctuations of collective variables. This can be done either by performing GCM 
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configuration mixing calculations of projected wave functions or by constructing 
collective Bohr-type Hamiltonians with deformation-dependent parameters deter- 
mined from self-consistent mean-field calculations. The possibility to perform self- 
consistent microscopic studies of shape transitions opens a new perspective on the 
origin of nuclear quantum phase transitions in various mass regions. It is there- 
fore important to systematically analyze these methods, also employing different 
energy-density functionals, various types of shape phase transitions that have been 
predicted in several regions of medium-heavy and heavy nuclei. 
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